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Abstract 



The electron electric dipole moment (dg) and the muon magnetic moment 
anomaly (a^) recently observed at BNL are analyzed within the framework 
of SUGRA models with CP violating phases at the GUT scale. It is seen 
analytically that even if de were zero, there can be a large Bino mass phase 
(ranging from to 27r) with a corresponding large B soft breaking mass phase 
(of size ~ 0.5 with sign fixed by the experimental sign of a^). The dependence 
of the B phase on the other SUSY parameters, gaugino mass mi/2, tan /3, 
is examined. The lower bound of determines the upper bound of mi/2- It 
is shown analytically how the existence of a non-zero Bino phase reduces this 
upper bound (which would correspondingly lower the SUSY mass spectra). 
The experimental upper bound on de determines the range of allowed phases, 
and the question of whether the current bound on de requires any fine tuning 
is investigated. At the electroweak scale, the phases have to be specified to 
within a few percent. At the GUT scale, however, the B phase requires fine 
tuning below the 1% level over parts of the parameter space for low and 
if the current experimental bound on de were reduced by only a factor of 3 — 4, 
fine tuning below 1% would occur at both the electroweak and GUT scale over 
large regions of the parameter space. All accelerator constraints (ruh > 114 
GeV, b — > S7, etc.) and relic density constraints with all stau-neutralino 
co-annihilation processes are included in the analysis. 
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I. INTRODUCTION 



The role that CP violation plays in elementary particle physics and how it relates to 
current theory, still remains after many years unclear. In the Standard Model (SM), CP 
violation is accommodated by inserting a single phase into the CKM matrix, and there are 
now under way a number of experiments to test the validity of this idea [0]. In supersym- 
metry, the CKM phase can also exist, but it is possible to have additional phases appearing 
in the soft breaking masses. From the viewpoint of string theory, CP violating phases are 
a natural occurrence. Thus in 10 or 11 dimensional M-theory models with six dimensions 
compactified on a Calabi-Yau manifold, the Kahler potential and Yukawa matrices are rep- 
resented by integrals over the Calabi-Yau space, and since this is a complex manifold, it is 
not surprising that CP violating phases can arise. With supersymmetry breaking, phases 
could arise also in the soft breaking masses. However, M-theory is not sufficiently developed 
to determine the details of such effects, and it is useful to have phenomenological constraints 
ClS db guide to where the fundamental theory may be. 

As is well known, the existance of CP violating phases in the SUSY soft breaking masses 
leads to electric dipole moments (EDMs) for the electron and neutron, and the smallness of 
the current experimental bounds on these puts severe constraints on the parameter space. 
One may satisfy these constraints by assuming that the phases are non-zero but small (i.e. 
0(10-2)) or that the SUSY masses are large (i.e. ~ 0(1 TeV)). The electron EDM (eEDM) 
arises from two diagrams, one involving the chargino-sneutrino intermediate state, and one 
involving the neutralino- selectron intermediate state. More recently it was suggested that 
a cancellation might occur between these two for relatively large phases (i.e. O(10~^)) over 
a reasonably large range of parameters |@], and there has been a large number of analyses 
based on this idea Diagrams similar to the above also enter into the muon g — 2, and 

the recently reported 2.6 a deviation from the Standard Model value for that quantity 
has placed significant bounds (at the 95% C.L.) on the allowed SUGRA parameter space [Q. 
It is thus natural to ask whether both constraints can be phenomenologically realized, and 
initial discussions of this have been made [^]. 

The minimal supersymmetric standard model (MSSM) has over 40 independent phases, 
and so cannot make significant theoretical predictions on this question. We use here instead 
supergravity (SUGRA) GUT models with gravity mediated supersymmetry breaking |T^ 
and R-parity invariance Q. Such models have become quite predictive, in part due to the fact 
that they automatically include the LEP results on grand unification, and because radiative 
breaking of SU{2) x U{1) implies that the SUSY soft breaking masses are of electroweak 
size. Thus the general size of both g — 2 and the EDMs (which depend on the size of the 
SUSY masses) are restricted. In addition, such models have a natural candidate for dark 
matter, the lightest neutralino (xi), and the condition that the relic density of neutralinos 
be in accord with the allowed range from astronomical measurements also puts important 



^Two alternate SUGRA models are anomaly mediated [^] and gauge mediated [12| soft breaking. 
The former appears to have difficulty in satisfying the Brookhaven E821 bounds on g — 2 when 
no SUSY CP phases are present jl^, while the latter does not appear to have a satisfactory dark 
matter candidate |14]. 
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constraints on the parameter space. 

Previous analyses of the EDMs within the framework of SUGRA GUT models has shown 
that a new fine tuning problem arises at the GUT scale for (the phase of B, the bilinear 
soft breaking mass) when tan/? gets large @]. Further, the discussion of the muon g — 2 
have shown that tan (3 is greater than 5 ( and very possibly greater than 10). In carrying 
out our analysis, then, we put a constraint on the parameter space that AObq/Obq > 0.01, 
where A^^o is the allowed range that will satisfy the experimental bounds on the EDMs. In 
order to carry out a complete analysis, however, it is necessary to include all the accelerator 
constraints (i. e. that rrih > 114 GeV (where h is the light Higgs boson), the CLEO bound 
on the 6 — > S7 decay, etc.) as well as the full analysis of the relic density including all the 
stau-neutralino co-annihilation channels. (The details of these were discussed in [|T5i[T6| .) 

Both the electron EDM and the muon g — 2 can be calculated from two types of diagrams, 
one with the intermediate chargino and sneutrinos states, and one with the intermediate neu- 
tralino and slepton states. If one assumes universal soft breaking in the first two generations, 
then the eEDM and the SUSY deviation of the muon magnetic moment from its Standard 
Model value, a^^^ {a^j, = {g^ — 2)/2), can both be calculated from a single amplitude for 
these diagrams, the former being related to the imaginary part, and the latter to the real 
part. In the following we will use the symbol to mean the deviation of the muon {g^ — 2)/2 
from its Standard Model value. 

As is well known, for large tan (3, is dominated by the chargino diagram |[T7|JT8|1 (the 
neutralino diagrams being generally quite small), while the strong experimental constraints 
on the eEDM require a near cancellation between the neutralino and chargino diagrams . 
In order to see how this can come about when CP violating phases are present, we calculate 
in the Appendix the leading large tan (3 terms. We discuss analytically in Sec. 2 how both 
experimental constraints can naturally be satisfied at the electroweak scale. In Sec. 3 we 
examine in detail numerically the combined experimental constraints of the eEDM and the 
Brookhaven E821 experiment on the SUSY parameter space for a variant of the mSUGRA 
model where the magnitudes of the soft breaking masses are universal at Mq but the phases 
are arbitrary. The experimental lower bound on puts an upper bound on the gaugino 
mass mi/2, and this upper bound is generally reduced when CP violating phases are present. 
The Higgs mass lower bound and the b ^ s'j constraint generally put a lower bound on 
mi/2- Imposing the fine tuning constraint at Mqut generally increases that bound when CP 
violating phases are present and also limits the range of the CP violating phases. Thus the 
two experiments interact to further restrict the SUSY parameter space when CP violating 
phases are large. In Sec. 4 we give some concluding remarks. 



II. MUON MAGNETIC MOMENT AND ELECTRON EDM 

We consider here supergravity models which are a generalization of mSUGRA to include 
the possibility of phases in the soft breaking masses at the GUT scale Mq = 2 x 10^^ GeV. 
The magnitude of the soft breaking masses are still assumed to be universal, but phases are 
not necessarily universal. Thus the SUSY parameters at Mq are mo (the scalar soft breaking 
mass), mi/2i = |mi/2| exp(z0j) i = 1,2,3 (the gaugino masses for the f/(l), SU{2), SU{3) 
groups), Aq = \Aq \ exp(iQ!o) (cubic soft breaking mass), Bq = \Bq \ exp^iOso) (quadratic soft 
breaking mass), and /io = l/^ol exp(i6'^) (the Higgs mixing parameter in the superpotential). 
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The model therefore depends on five magnitudes and six phases. However, one can always 
set one of the gaugino phases to zero, and we chose here 02 = 0. 

The renormalization group equations (RGEs) allow one to evaluate the parameters at 
the electroweak scale. To one loop order, the 0, and 9^ do not run. Further, radiative 
breaking of SU{2) x f/(l) at the electroweak scale allows one to eliminate 6^ in terms of 
6b, the B phase at the electroweak scale, and determine and \B\ in terms of the other 
parameters and tan/5 = \ {H2) / {Hi)\. Thus with a convenient choice of Higgs VEV phases, 
the minimization of the Higgs potential yields ||19|| : 



0^ = -Ob + f\-eB + -Ob + ai) (1) 
, ,9 m? — tan^/3mn , ^, 1 ^ ml 

where ag and are the quark and lepton phase of Ag and Ai at the electroweak scale, 
mf = m?jj_ + (z = 1,2), m| = 2|yUp + m\ + m^, and m^/- are the Higgs running masses 
at the electroweak scale. Sj = dVi/dvf where Vi is the one loop contribution to the Higgs 
potential, and Vi = \ (Hi) \ . There remain therefore four real parameters which we take to be 
mo, l^ol) 1^^1/21 and tan/5, and four phases 6*^05 cto, (pi and 03. In this paper we examine 
only the electron electric dipole moment^ (eEDM) , and so our results are only very weakly 
dependent on 03. Also, since we are dealing only with first and second generation sleptons, 
there is only a weak dependence on oq. Thus the two important phases are 6bo and 0i. The 
parameter range that we examine is: 

|mo| < 1 TeV; \my2\ < 1 TeV; \Ao/mi/2\ < 4 (3) 



tan /3 < 40 (4) 

As discussed above, both and dg can be obtained from the same complex amplitude 
A (assuming universal soft breaking in the first two generations). One has then 

■mlRelA] (5) 



dele = ^— — melmMl (6) 

' STTsin^^vF " ^ J ^ ^ 

The amplitude A is defined in the Appendix. For the case where there are no CP violating 
phases, the experimental data favors large tan /3 |^ . In order to see semi-quantitatively 
the nature of the cancellations needed to make de small, the leading terms for large tan/? 
were obtained in the Appendix when /i^ ^ (which is almost always the case for the 



^While there has been considerable progress in calculating the neutron EDM 
EDM there remain still significant hadronic uncertainties, in contrast to the clear calculation 
of the eEDM. Further, the 03 phase can always be adjusted to satisfy the neutron EDM, and the 
effect of the ^^^Hg EDM would then only result in reducing the remaining parameter space that 
we find from the eEDM and the muon g — 2 given here. 



2C| and ^^^Hg 
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mSUGRA parameter space). Prom Eqs. (A18, A25) we find for the neutralino and chargino 
diagrams 



where 
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a m^ tan /3 
47rm2|A*| sin^ 



(10) 



and C^(x°) = tan^ ^^vF(^2/|^i|)CAt(x^)- The form factors Fj = F(m?/m?±) and G^j = 
(^(m^^/mlo) are defined in Eqs. (A8, A12), and |mi| = 0.4mi/2, ''tT'2 — 0.8mi/2. (Our states 
are labeled such that e.g. m^o < m^q for i < j.) 



A similar decomposition of the electron electric dipole moment, dg/e = d{x 



gives: 
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where = ('^e/2m^)C^. 

In the case where all phases are zero, it is well known that the chargino contribution to 
dominates over the neutralino diagram even though the front factors C^(x^) and C^(x°) 

are of comparable size. This can be understood from Eqs. (8) and (9) in the following way. 

The second term in Eq. (8), coming from the heavy chargino, xt^ cancels about 30% of the 

light chargino contribution. In a similar fashion, the second term in Eq. (10) coming from 
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the heavy smuon , jj,2, contribution cancels part of the leading term. However, due to the 
slowly varying nature of the form factors Gu and G21, about 75% of the leading Gu term 
is canceled, reducing a(x°) significantly. The remaining two terms arising from the heavy 
neutralinos, X234 ^'^^ generally small (and there can be cancellations also between these 
terms). In contrast, when the CP violating phases are not zero and an electric dipole moment 
exits, the neutralino and chargino contributions must be of nearly equal size and opposite 
sign so that dg be greatly suppressed to be in accord with experiment. This change in the 
relative sizes of the two terms must be brought about by the sine factors in Eqs. (11, 12). 
Simultaneously, the corresponding cosine factors in Eqs. (8, 9) will modify the predictions 
of a^. This then leads to two questions: 

(1) Can the experimental constraints on dg and (and of course all other experimental 
constraints) be simultaneously satisfied with "large" phases, i. e. of O(10~^)? 

(2) If large phases are possible (and d^ is very small) can the experimental constraints 
be satisfied without undue fine tuning of the phases? 

Question (2) divides into two parts: (2a) Is there undue fine tuning of phases needed to 
satisfy the experimental constraints at the electroweak scale, and (2b) Is there undue fine 
tuning of the phases needed at the GUT scale? In the following, we define the fine tuning 
parameter R{(j)) for any phase by 

= 1/ , , , X = -— (13) 

2(^2 + 01 ) 0av 

where 02 and 0i are the upper and lower value of when the experimental constraints on 



de and are both satisfied. We use here the current bound on dg |]21| : 

|4| < 4.3 X 10"^^ e cm (14) 

and the 2 std range for 0: 

11 X 10"^° < < 75 X 10"^° (15) 
In the following we will assume that no fine tuning less than 1% be allowed, i. e. 

i?(0) > 0.01 (16) 

Within this framework, we find that Question (1) can be answered affirmatively, i. e. 
cancellations between the neutralino and chargino diagrams in de can indeed be satisfied 
with large phases. Further, the answer to Question 2(a) is that no undue fine tuning of 
the electroweak parameters is necessary to achieve the suppression of d^. However, we will 
see that this is not the case at the GUT scale, where the simultaneous requirements of 
electroweak radiative breaking and the experimental d^ constraint leads to significant fine 
tuning, and if Eq. (16) were imposed, a large amount of the parameter space is eliminated. 

In a GUT theory, the fundamental parameters are specified at Mqut? and the conse- 
quences at the electroweak scale are obtained from the renormalization group equations 
(RGEs). These GUT parameters are presumably to be determined at some future time by 
a more fundamental theory (e. g. string theory), and so fine tuning at the GUT scale may 
imply a significant theoretical problem (while fine tuning at the electroweak scale would be 
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an acceptable theoretical consequence of the RGEs). Of course, what level of fine tuning 
is acceptable is somewhat a matter of choice, and we view Eq. (16) only as a reasonable 
benchmark to consider. 

In Sec. 3 below, we will consider these results in detail quantitatively. We give here an 
analytic discussion of how they arise. To show that large phases are easily achievable, we 
write de in the form 

4/e = -D,{x^)A[siYi{e^) + asin(^^ + 00] (17) 

where the coefficients A and a can be read off from Eqs. (11, 12). In the extreme case where 
4 = 0, Eq. (17) and Eq. (1) imply 

asin</)i 

tanfc'B = (18) 

1 + a cos 01 

where we have neglected the small 1-loop correction in Eq. (1). The fact that the chargino 
diagram dominates over the neutralino diagram implies that a < 1, and detailed numerical 
calculations show that a ~ 0.2 — 0.4 for much of the SUSY parameter space. Thus as 0i 
varies from to 27r, over most of the parameter space \9b\ will be large (rising to a maximum 
of about 0.5) even though d^. has been set to zero! 

We next consider the effects of the CP violating phases on a^. From Eqs. (8, 9), we have 

a^, = C^{x^)A[cos{e^) + a cos{e^ + 0i) + b] (19) 

where b is the term in Eq. (9) independent of the phases. In general, b is quite small, and 
we will neglect it in the following. (Of course, in the numerical calculations in Sec. 3 all 
such effects as well as the loop corrections are considered.) Using Eq. (18) to eliminate 9b, 
we obtain 

= ±a^(0) g(a,0i) (20) 
where a^(0) is the value of with zero phases, 

^~ (1 + a) ^^^^ 

The ± factor in Eq. (20) is the sign of cos6b- Since experimentally, > 0, this implies 

eB>0 for < 01 < vr; Ob <0 for vr < 0i < 27r (22) 

The two branches of Eq. (22) are symmetric, and in the following we consider the > 
branch. The factor Q in Eq. (21) reduces the theoretically expected size for a^. Since the 
experimental lower bound on implies an upper bound on mi/2 [§], the Q factor due to 
the phases will reduce this upper bound, further restricting the allowed SUSY parameter 
space. However, since in general Q ~ 0.5, this reduction will still be consistent with all 
experimental data, and the effect of the SUSY CP violating phases does not qualitatively 
change the fit to the data for that was obtained assuming no CP violating phases. 

We can also estimate how much fine tuning is needed in the phases to satisfy the exper- 
imental bound of Eq. (14). Thus let A6b be the change in 6b for a fixed value of 0i as de/e 
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varies from —4.3 x 10 cm to +4.3 x 10 cm. Characteristically, the factor D{x^)A in 
Eq. (17) is numerically about 100 times the current upper bound on df. (e.g. for mi/2 = 480 
GeV (mo = 118 GeV), = 690 GeV and tan/3 = 15, this factor is 4.1 x lO'^^ cm.) 
Considering then the variation of in Eq. (17) one finds 

2 X 10"^ ^ A9b cos{9b) [1 + a cos(0i) + a sin((/)i) tan(^B)] (23) 

Hence using Eq. (18), 

R{eB) = AOb/ sm{eB) = 2 X 10"V(« sin(0i)) (24) 

Thus the allowed range of 6b is indeed small, though the condition of Eq. (16) is generally 
satisfied. A reduction of d^. by a factor of 10 would be enough to produce a serious fine 
tuning problem at the electroweak scale. In contrast, we will see below that there is already 
a significant fine tuning problem at the GUT scale even with the current bound on de- 



III. DETAILED CALCULATIONS 

In this section, we consider detailed calculations of the effects of the experimental con- 
straints involving the eEDM and which were analytically estimated in Sec. 2. The analysis 
is done within the framework of the generalized mSUGRA described in Sec. 2 (though exten- 
sion to non- universal soft breaking models can easily be done). In order to get a clear picture 
of what constraints on the SUSY parameter space arise, it is necessary to simultaneously 
impose all the experimental constraints. Aside from the above, these include the following: 
(1) The LEP Higgs mass lower bound rrih > 114 GeV [^. Since the theoretical analysis of 
mh still has about a 3 GeV uncertainty, (which (conservatively) may be an overestimate) 
we interpret this in the theoretical calculation |TH] to mean rrih > 111 GeV. (2) The b ^ sj 



branching ratio constraint. We take here a 2 std range of the experimental CLEO data |23 



1.8 X 10"^ < BR{B Xs7) < 4.5 x 10"^ (25) 
In the theoretical analysis we include the NLO SUSY contribution for large tan/? [^, as 



these produce significant effects (particularly since the lower bound favors larger values 



of tan/5). (4) We include the 1-loop corrections to the b and r masses ||2^, which are 
significant also for large tan j3. (5) All co-annihilation effects are included in the relic density 
calculations. We assume here the range for the neutralino cold dark matter to be 

0.025 < n^oh^ < 0.25 (26) 

The upper bound is consistent with recent Boomerang and Maxima measurements |]26| , |27| , 
while the lower bound allows for the possibility that there may be more than one type of 
dark matter. However, our calculations here are insensitive to the precise value of the lower 
bound, and one would get very similar results if the lower bound were raised to 0.05 or 0.1. 

The effects of the Higgs mass and b —>■ S'j bounds is to push the allowed parameter 
space mostly into the co-annihilation domain of the relic density calculation. (Thus only 
a small part of the allowed parameter space occurs at small enough mi/2 to lie below the 
region of co-annihilation, which begins at mi/2 — 350 — 400 GeV.) In SUGRA models, stau - 
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neutralino co-annihilation can occur quite naturally. Thus for mo = 0, the charged sleptons 
lie below the lightest neutralino, Xi; and one must increase mo to raise their masses so that 
the Xi is the dark matter particle. Thus as mi/2 increases, tjiq must be increased in lock step 
so that the neutralino remains the LSP, and one finds a relatively narrow corridor in the 
mo — mi/2 plane consistent with the relic density bound of Eq. (26) and with the lightest stau 
lying above the neutralino. The dependence of these corridors on tan/5 and Aq are shown 
in Figs. 1 and 2 for the case where there are no SUSY CP violating phases. We see that the 
corridors (which are characteristically 20-30 GeV wide) lie higher for larger tan (3 and larger 
|y4o|. This is because the stau mass decreases when these parameters are increased, and so 
one must raise mo to keep the stau mass greater than the neutralino^. In general, the effect 
of co-annihilation is to determine mo approximately in terms of mi/2 for a given tan/3 and 
Ao. This greatly sharpens the predictions of the theory. 
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400 
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300 400 500 600 700 800 900 1000 
FIG. 1. Corridors in the mo — mi/2 plane allowed by the relic density constraint for (bottom 
to top) tan/? = 10, 30, 40 for m/j > 114 GeV, Aq = 0, fi > 0, all phases set to zero. The tan/3 = 30 
and 40 corridors all lie in the co-annihilation region, while only the beginning of the tan /? = 10 
corridor is in the non co-annihilation region. The Higgs mass bound determines the lower mi/2 
bound for tan /3 = 10, while both the Higgs mass and the b sj bounds equally produce the lower 
bound for tan/3 = 30 and b ^ determines the lower bound for tan/3 =40. The short slanted 
lines represent the upper bound on mi/2 due to the lower bound, of Eq. (15). 



A. Allowed Regions at the Electroweak Scale for 6b and 01 

We now examine the effects of having non-zero CP violating phases present, and discuss 
the dependence of the allowed phases on the SUSY parameters. To illustrate the phenomena, 
we consider one low tan f3 and one high tan /?. Since for no CP violating phases, one had (at 



^Increasing tan/3 or making ^o negative increases the magnitude of the L-R term in the stau 
mass matrix of Eq. (A3) and hence decreases the light stau mass. For Aq > 0, the opposite effect 
occurs, but also the diagonal matrix elements of the mass matrix are reduced, and again the stau 
mass decreases. 
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400 500 600 700 800 900 1000 

mj (GeV) 

2 

FIG. 2. Corridors in the mo — mi/2 plane allowed by the relic density constraint for (from 
bottom to top) ^0 = 0, — 2mi/2) 4mi/25 for rrih > 114 GeV, ^ > 0, all phases set to zero. The 
lower bounds for Aq = 0, —2mi/2 are due to the b ^ constraint, and for Aq = 4mx/2 from 
the Higgs mass bound [15]. 

90% CL.) tan/5 > 10, we examine the cases of tan/5 = 15 and tan/5 = 40. Fig. 3 shows the 
corridors allowed for 6b (the B phase at the electroweak scale) by the current experimental 
bounds on de for tan/5 = 40, Aq = for two choices of the gaugino phase: 0i = 0.9 (lower 
curves) and 0i = 1.2 (upper curves). One sees that one gets a significant phase 9b of the 
size expected by Eq. (18), the larger 0i allowing a larger 9b, also in accord with Eq. (18). 
Note also that the allowed corridors for 9b widens as mi/2 increases as expected, since the 
experimental constraint is less severe for a heavier mass spectrum. A similar graph is 
shown in Fig. 4 for 0i = 0.9 (upper curves), 0i = 3.4 (lower curves). 9b turns negative for 
01 in the third quadrant, again as expected from Eq. (22). 

We note in all these curves, the upper bound on mi/2 (due to the lower bound on a^) 
is reduced compared to the case when the CP violating phases are zero. (Then the upper 
bound is mi/2 = 790 GeV for tan /5 = 40, Aq = This is due to the phase 0i in the Q 
factor in Eq. (20). Q is smallest when (pi is near vr, as can be seen in explicitly in Fig. 4 for 
</>! = 3.4 

We consider next the tan/5 dependence of the allowed region for 9b- This arises due 
to an indirect tan/5 dependence in the parameter a of Eqs. (17,18). As seen from Fig. 2, 
coannihilation determines mo in terms of mi/2, and this mg increases with tan/5, changing 
the ratio of neutralino to chargino diagram differentially. Fig. 5 shows the allowed region for 
tan/3 = 15, and Aq = 0, 0i = 1.2 (upper curves) and 0i = 0.9 (lower curves). We see that 
one can get considerably larger values of 9b at lower tan /?, though the upper bound on mi/2 
due to the lower bound on is considerably reduced at low tan /?. The Aq dependence is 
exhibited in Fig. 6, where the allowed corridor for 9b is plotted for tan/5 = 40, (pi = 0.9 and 
Aq = (upper curves), |y4o| = 2mi/2, = 0.5 (lower curves). Increasing the magnitude of 
Aq increases the value of uiq (by Fig. 2) and reduces the size of 9b- 
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FIG. 3. Regions allowed for Ob hj the experimental constraint on as a function of foi' 
tan/? = 40, Aq = 0, for cpi = 0.9 (lower curves) and 0i = 1.2 (upper curves). 
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FIG. 4. Same as Fig. 3 for (pi = 0.9 (upper curves) and (pi = 3.4 (lower curves). 

B. Fine Tuning at the Electroweak Scale 

The above analysis shows that the dg experimental constraint allows 6b to be 0(10^^) 
for a wide region of parameter space, and 0i can range widely, i.e. from to 2tt. We 
next investigate whether the smallness of the upper bound on d^. requires any fine tuning to 
maintain this constraint. In Fig. 7 we plot the fine tuning parameter R of Eq. (13) for 6b 
for tan/? = 40, Aq = for (pi = 0.9 (upper curve) and (pi = 1.2 (lower curve). One sees that 
R{6b) is small, i.e. a few percent, but satisfies the criteria R > 0.01 for the entire range. The 
size of R{6b) is consistent with what was expected from the analytic analysis of Eq. (24), 
and increases as (pi decreases also as expected. Fig. 8 shows a similar plot for tan/? = 15. 
Since 6b is larger here (see Fig. 5), R{6b) is smaller, but still within the acceptable range. 
Fig. 9 shows R{(pi) for tan/? = 40, Aq = for 6b = 0.2 (upper curve), and 6b = 0.3 (lower 
curve). Again R{(pi) > 0.01 for the entire range of mi/2, and is generally larger than R{6b) 
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FIG. 5. Same as Fig. 3 for tan/? = 15. The allowed regions terminate at low due to the 
m/j constraint. 
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FIG. 6. Allowed region for 9b for tan/5 = 40 for = (upper curves) and |ylo| = 2mi/2) 
CKo = 0-5 (lower curves). 

(by about a factor of tan(/)i/(/)i). 

We see from the above results that the current experimental bound on df, can be accom- 
modated with large phases. Further, for most of the parameter space, while the range of the 
allowed phases at the electroweak scale required to accommodate the bound on de is small, 
no fine tuning below 1% is needed. We will see, however, a more serious fine tuning problem 
develops at the GUT scale. 

C. Fine Tuning at the GUT Scale 

We now examine what parameters are fine tuned at the GUT scale to achieve the experi- 
mental EDM bounds. From Eq. (1), we see since the loop correction is small, that 6'^ = —6b, 
and since we have seen that 6b does not need excessive fine tuning, the same can be said for 
6^. Further, since 6^ does not run with the RGE, its value at the GUT scale is the same as 
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FIG. 7. RiOs) as a function of fo'^ tan/? = 40, Aq = 0, </>i = 0.9 (upper curve), 0i = 1.2 
(lower curve). 
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FIG. 8. Same as Fig. 7 for tan (3 = 15. 

at the electroweak scale and hence no fine tuning is needed at Mq- A similar result holds for 
01, which also does not run with the RGE. However, as has been previously pointed out [Q, 
matters are different for the B phase at M^, 6'^^, and we review briefly the discussion given 
thereQ. To see analytically what is occurring, we consider the intermediate and low tan/? 
region, where the RGE can be analytically solved. One flnds for B the result 

5 = fio - ^(1 - ^o) -Y.^^\my^\e'^^ (27) 

where Dq = 1 — imtj sin/5)^ ~ 0.25 and $j = 0(1). As tan/3 gets large, the radiative 
breaking condition Eq. (2) shows that \B\ gets small. Taking the imaginary part of Eq. (27) 
gives 





^This analysis differs from that given in |25] which does not take into account the possibility of 
cancellations in between the neutralino and chargino diagrams. In fact the discussion in l2J 



sets the phase to zero. 
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FIG. 9. R{4>i) as a function of for tan/3 = 40, Aq = Q and 6b = 0.2 (upper curve), 
6b = 0.3 (lower curve). 

|5|sin^B = l^olsin^B, - (1/2)(1 - L>o)|^o| sinao - I] <^'^|"^l/2| sin0i (28) 

To the lowest approximation, since |-B| is small, one may then neglect the Ihs of Eq. (28). 
Eq. (28) may then be viewed as an equation determining Obq in terms of the other GUT 
scale phases, and hence 9b(^ will in general be large if the other phases are not all small (a 
result that is confirmed in by detailed calculation). However, the range of 9bq so that the 
experimental bound on de is satisfied is then significantly reduced. Thus for fixed values of 
ao and 0j, Eq. (28) gives as tan/3 gets large (i.e. \B\ becomes small) 

MB, = mi\B^\)^OB<:MB (29) 

Since we have already seen that AO b is small (though not violating the fine tuning condition) , 
one may expect that AObq may need to be fine tuned, particularly for low mi/2 where the 
SUSY mass spectrum is light. This is seen explicitly in Fig. 10 for tan/5 = 40, = with 
(pi = 0.9, 1.2, 1.6, 2.3, and 2.6. We see that R{Obo) decreases as 0i increases from 0.9 to 1.6 
(= 7r/2) and then increases for 2.3 and 2.6 where tt — 01 is decreasing. (The upper bound 
on mi/2 arising from the lower bound on decreases as (pi moves into the second quadrant 
in accord with Eqs. (20,21).) We see that if one imposes the fine tuning constraint that 
R{6bo) > 0.01, large sections of the low mi/2 region would be excluded, e. g. one would 
require mi/2 > 540 GeV for (pi = 1.6. The fine tuning becomes more acute at lower values of 
tan/3. Thus Fig. 11 shows R{Obq) as a function of mi/2 for tan/3 = 15, = for (pi = 0.9 
(upper curve) and (pi = 1.2 (lower curve). Thus the fine tuning constraint would eliminate 
completely (pi = 1.2 (and the entire region of ~ 0.4 radians around vr/2) and also restricts 
the other values of (pi. Finally, we note that increasing Aq generally increases the amount 
of fine tuning needed. Thus Fig. 12 shows R{6bq) for tan/? = 40, (pi = 0.9 for Aq = 
(upper curve) and |y4o| = 2mi/2, Oo = 0.5 (lower curve). The entire |y4o| = 2mi/2 curve has 
R{9bo) < 0.001. 
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FIG. 10. R{0Bo) as a function of mi/2 for tan/3 = 40, = for (from bottom to top) 
01 = 1.6, 1.2, 0.9, 2.3 and 2.6. 
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FIG. 11. R{9bq) as a function of mi/2 for tan/3 = 15, tIq = for (pi = 0.9 (upper curve) and 
= 1.2 (lower curve) 

IV. CONCLUSIONS 

We have examined here what SUSY CP violating phases are possible within the frame- 
work of SUGRA models, when the electron electric dipole moment experimental bound is 
imposed on the SUSY parameter space. In order to analyze this, we imposed simultaneously 
the accelerator bounds of the Higgs mass {mh > 114 GeV) and the b —>■ s'y branching ratio, 
the relic density bound for neutralino cold dark matter and the recent 2.6 std deviation of 
the muon magnetic moment from the Standard Model prediction. The different experimen- 
tal constraints interact with each other. Thus the Higgs mass and 6 — 57 constraints put 
lower bounds on the gaugino mass (mi/2 ~ (300 — 400) GeV) which puts the relic density 
analysis mostly in the region where stau-neutralino co-annihilation occurs. This closely fixes 
the scalar mass ttiq in terms of mi/2 (for fixed Aq and tan/?). The lower bound then puts 
an upper bound on mi/2- Thus the parameter space becomes highly constrained. One can 
estimate analytically, as was done in Sec. 2, the effects of turning on the CP violating 
phases. In fact if were zero, one could still have large CP violating phases present, with 
the Bino phase (pi between and 27r. The condition that be positive puts the B phase 
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FIG. 12. R{6bq) as a function of for tan/? = 40, (pi = 0.9 for Aq = (upper curve) and 
Ao\ = 2mi/2; tto = 0-5 (lower curve). 



9b in the first and fourth quadrants with \6b\ ~ 0.2 — 0.4. The dependence of 6b on tan/? 
and was discussed in detail in Sec. 3. It was seen that the phase (pi acts to reduce the 
theoretical value of by a factor Q < 1, defined in Eq. (20, 21), and from the experimental 
lower bound on this reduces the upper bound on mi/2, limiting further the parameter 
space. 

The relevant question, however, is whether the smallness of the experimental bound 
on de requires an unreasonable amount of fine tuning of the phases. Using the parameter 
R{(j)) = A0/0av, we find at the electroweak scale, both R{6b) and R{4>i) are small, i. e. 
a few percent. However, in general R > 0.01, and so no fine tuning below the 1% level is 
needed. However, at the GUT scale, this is not the case for R^Ob^) for a significant part 
of the parameter space. Thus if we were to exclude regions where R{6bq) < 0.01, then for 
tan/5 = 15, Aq = 0, (pi phases near 90° ( i. e. 1.2 < 01 < 2.0) are completely excluded 
(Fig. 11). The effect is reduced for higher tan/5. However, for example, for tan/5 = 40, 
Aq = the condition R > 0.01 would eliminate mi/2 < 540 GeV for (pi = 1.6, and raise 
the lower bound on mi/2 by a lesser amount for other values of (pi (Fig. 10). Increasing Aq 
decreases the value of R, making the fine tuning more serious, as can be seen in Fig. 12. 

The fine tuning problem is thus on the verge of becoming quite acute. The experimental 
bound on de is likely to decrease by a factor of three in the near future |2^. The effect 
this would have is shown in Figs. 13 and 14. Fig. 13 shows R{6b) as a function of mi/2 for 
tan/5 = 40, Aq = 0, 0i = 0.4 (upper curve) (corresponding to 9b — 0.1) and 0i = 0.9 (lower 
curve). The curves are what would occur if the experimental bound on de were reduced 
by a factor of three. Fig. 14 shows R{9bo) as a function of mi/2 for tan/5 = 15, Aq = 0, 
(pi = 0.3 (upper curve) (corresponding to 9b — 0.1) and (pi = 0.9 (lower curve), if the bound 
is reduced by a factor of three, (pi < 0.4 in Fig. 13 and (pi < 0.3 in Fig. 14 would have less 
fine tuning but correspond to 9b < 0.1. 

In a GUT theory, presumably the parameters at the GUT scale are the more fundamental 
ones, and fine tuning of these parameters would presumably represent a serious problem. 
Of course, what level of fine tuning one accepts is somewhat a matter of taste, and we 
view the criteria R > 0.01 as merely a bench mark for consideration. However, the above 
results show that any significant experimental reduction of de would make the idea of large 
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FIG. 13. R{Ob) as a function of mi/2 for tan/3 = 40, Aq = 0, = 0.9 (lower curve) and 
= 0.4 (upper curve)with the current experimental bound on de reduced by a factor of three. 
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FIG. 14. R{9bq) as a function of mi/2 for tan/? = 15, Aq = 0, (/>i = 0.9 (lower curve) and 
(f)i = 0.3 (upper curve)with the current experimental bound on de reduced by a factor of three. 

SUSY CP violating phases more difficult to maintain, as fine tuning would set in even at 
the electroweak scale unless the SUSY masses are heavy. 

This work was supported in part by National Science Foundation grant No. PHY- 
0070964. 

APPENDIX A: DERIVATIONS 

In this Appendix, we calculate the leading terms in tan (3 for the complex amplitude A 
of Eqs. (5) and (6) needed to calculate and d^. In order to do this, it is necessary to 
diagonalize the various mass matrices entering in the chargino and neutralino loops. These 
matrices depend on the phases 0i, ai and 9 = 9^ + ei + €2 where ai is the phase of Ai and 
(ifi_2) = ^^1,2 e**^^'^ where Vi^2 = |(-f^i,2)|- Minimizing the effective potential with respect to 
ei^2 determines 9 in terms of 9b, and one may then chose Higgs phases such that ei 2 = 0, as 
was done in Eq. (1). With this choice of phases, the chargino and neutralino mass matrices 
are 
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m2 -\/2Mvi/ sin (3 
V^Miy cos (3 /i 



(Al) 



/ mi a 6 

m2 c (i 

a c — /i 

V & d -/i 



(A2) 



where fhi = |mi|e"'^i, /i = |/i|e a = 
d = —6 cot 6'vi/, and m2 has been chosen real and positive. 
The slepton mass matrices have the form 



Mz sin 9w cos P, b = Mz sin 9w sin /3, c 



-a cot ^ly. 



ILL 
,2 * 



(A3) 



where mf^,^ = m;(y4f — /itan/5), A; 



m 



ILL 



2 I 2 

rrii^ + nil 



\Ai\e 
1 



and 



1 — 2 sin ^w')M|cos 2(3 



m 



IRR 



m\ + mf — sin^ O^^M"^ cos 2/3 



(A4) 
(A5) 



The m| ^ are obtained by running the RGE from the GUT scale to the electroweak scale |pO|| . 
In our analysis we will assume universal soft breaking of the ft and e scalar masses (mg) and 
universal cubic soft breaking masses (Aq) at Mq. Since ml and are very small, they can 
be neglected at the electroweak scale, i.e., M? = M|. 

M^o is a symmetric, complex matrix and can be diagonalized by a unitary matrix X 
according to M^oX = X*M^^^ where 



diag(m^o, m^o, m^o, m^o] 



(A6) 



Mj can be diagonalized by a hermitian matrix D with M^D 



DAf^^^^ where M~ 



diag(m? 



2(D) 



Finally one diagonalizes M^± by two unitary transformations U and 



M^±^ where M^f ^ 



diag(m-±,m-±). 



according to U*M^±V^ 

The amplitude A can be divided into its chargino and neutralino parts: A = A{x^] 
A{x^). We follow the notation of PT| where one finds that 



+ 



A{x 



1 



V2Mw cos (3 



(A7) 



and F,- 



F(m?/m2±; 



with 



Fix) 



1 — 3a; 



X) 



2x^ In X 

(1-X)3 



(A8) 



Similarly 
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(A9) 



where 



n1 



—j={taii9wXij + X2j)D\^ — KiX3jD2i. 



X 



(AlO) 
(All) 



V2tanewXijD2k + niXsjDik 
X!; = ^tan^ew\Xij\^[\Dik\^ + A\D2k\^) + ^\X2j\^\D 
+ tan 9w \Dik\ ^XijX^j + 0{mi) 

and Ki = mi/ {\/2Myy cos/5). The loop integrals are G^j = G{mf /m|o), H^j = H(mf /m\o) 
where 

(A12) 



G(x) 



1 + X 2x 

+ — 7^ Inx 



n-x] 



(1 — x)-^ 



(1 -X)3 

6x 
(l-x) 



■ \nx 



(A13) 



We need only keep the terms in A independent of mi and thus can neglect the 0{mi) terms 
in Eq. (All). (Actually, rif begins linearly in mi.) 

We are interested in calculating only the leading terms in tan (3. We will also do this 
in the limit M^/|/ip <^ 1 and M^/|mjp <^ 1 (which is valid for most of the mSUGRA 
parameter space). In that case one has [0] 



m^o = mi : m^o = m-± = mo; m^o = m-± = u 



(AM) 



We consider first the calculation of Alx"^)- Since V diagonalizes M^±M^± and U* diagonal- 



izes M^±M^± one finds for the leading terms 



TT* — 
^12 — 





1/^ 


2 


l/^l 


2 _ 


- ml 



V21 = -^V2Mwsm/3--P^V, 



(A15) 

. , .22 (A16) 

m2 — mg 

With an appropriate choices of phases one has to lowest order Uu = 1 = U22, Vu = 1 and 

V^22 = e''^ (A17) 

Hence inserting into Eq. (A7) gives 

tan P 



m2|/i| 





l/i 


2 




2 _ 


- ml 



m^ 

I 10 ~ 9 ' 

l/ip — m2 



(A18) 



To calculate the leading terms of A{x^) it is useful to first note the size of the matrix 
elements Xij. Thus to zeroth order in Mz 
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X22 — 1 



L44 



e 2 X 



33 



(A19) 

(A20) 



Also one has X12, X21 are 0(M|) (and hence neghgible) while the remaining components 
are 0{Mz). To lowest order, the slcpton mass eigenvalues are = 'fn'iRR-i f^j^ — ''^'ill-i ^^<i 
the D matrix has the form D12 = 1 = D21 and 



mi{A*i — //tan/3) 

2 2 
mf — mf 



-D 



22 



(A21) 



To illustrate the calculation of A{^) we consider the leading term when k = 1 
Prom Eq. (AlO), two terms contribute to rjl for large tan/3: 

Vl = - (^^tan^Xni^tij (v^tan^H^Xn£>2i) 
+ (/^iXaiD^i) (V2tan^w'Xni?2i) 



which evaluates to 



1 mi tan^ 6w tan /3 
?7i = 



where we have used 



^31 = 



2 2 

— 

12 (1 





/i 


2 


1/^ 


2 




mi 


2 



(A22) 



(A23) 



Mz sin 6*^^ sin /? 




2 






2 




mi 


2 
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(A24) 



(Note that ril is linear in rrii and hence Eq. (A9) is not singular as mi 0.) In a similar 
fashion one can obtain all the leading terms in Eq. (A9). (The Xj terms of Eq. (All) do 
not contribute.) The total answer is 

tan^ 9w tan P 



Am 



|mi||/i| 

2 2 
mf — mf 

<2 tl 



2 2 

mf — mf 

h 11 





^1 


2 


l/i 


2 




mi 


2 
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1)" 


2 




mi 


2 



1 1 
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2 tan^ 9w 




2 _ 


- mi 



mi 



G22 — - 



^21 je 
1 |mi|m2 



G23 e-*^" 



1 1 1 


fhi 






2 


4 tan^ 6w 








2 _ 


- ml 



(A25) 



We note that a large amount of cancellation occurs in this regime: terms proportional to 
G24 have all canceled with part of the G23 terms, and the total G14 contribution cancels with 
the Gis terms. Note also that the A{^) depends separately on two phase combinations: 
9n + 4>i and 9^, though terms depending on 9^ + 4>i are generally larger. 
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